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At  the  end  of  the  first  century  BC,  Ephesus  became  the  Roman  capital  of  Asia  Minor  and  the  most 
important  commercial,  religious,  and  cultural  center  of  the  region.  In  order  to  evaluate  the  status  of 
anthropogenic  fluxes  in  the  port  of  Ephesus,  a  12  m  long  sediment  core  drilled  in  the  Roman  basin  was 
investigated  to  shed  light  on  the  paleo-environmental  evolution  of  the  harbor  using  grain  size  distri¬ 
bution  analysis,  14C  ages,  major  and  trace  element  geochemistry,  and  Pb  isotope  compositions.  With  the 
help  of  complementary  sedimentological  data  and  Principal  Component  Analysis,  five  distinct  units  were 
identified  which,  together,  reflect  the  different  stages  of  water  history  in  the  harbor.  Among  the  major 
disruptive  events  affecting  the  port  were  earthquakes  and  military  events,  both  of  which  were  partic¬ 
ularly  effective  at  destroying  the  water  distribution  system. 

Seasonal  floods  of  the  Cayster  River  (Kugiik  Menderes)  were  the  major  source  of  the  silt  that  progres¬ 
sively  infilled  the  harbor.  Silting  in  was  further  enhanced  by  the  westward  migration  of  the  river  mouth.  A 
single  major  disruptive  event  located  at  550  cm  core  depth  and  heralding  the  development  of  anoxia  in  the 
harbor  marks  the  end  of  the  dynamic  regime  that  otherwise  controlled  the  harbor  water  throughout  the 
Roman  Empire  period.  This  remarkable  event  may  correspond  to  a  major  disruption  of  the  aqueduct  system 
or  to  a  brutal  avulsion  of  the  Cayster  River  bed.  It  clearly  represents  a  major  disturbance  in  the  history  of  life 
at  Ephesus.  It  is  poorly  dated,  but  probably  occurred  during  the  reign  of  Augustus  or  shortly  after.  Lead 
isotope  and  trace  metal  evidence  suggest  that  in  the  four  bottom  units  pollution  was  subdued  with  respect 
to  other  Pb  metal  inputs,  presumably  those  from  aqueducts  and  natural  karstic  springs.  Near  the  top  of  the 
core,  which  coincides  with  harbor  abandonment  and  the  more  recent  period,  anthropogenic  Pb  contam¬ 
ination  is  clearly  visible  in  both  Pb  abundances  and  isotopic  compositions. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Lead  isotope  studies  have  opened  up  a  new,  though  somewhat 
controversial,  perspective  on  the  development  of  the 
manufacturing  status  of  ancient  cultures  over  the  past  several 
millennia  (Hong  et  al.,  1994).  Isolated  artifacts  alone  do  not  suffice 
to  assess  the  broad  and  long-lasting  aspects  of  antique  trade  routes. 
Lead  isotopes  constitute  a  complementary  tool  in  that  they  play  a 
critical  role  wherever  their  compositions  can  be  ascribed  to 
anthropogenic  influence  in  the  form  of  lead  and  heavy  metal 
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pollution  of  sediments  accumulated  in  harbors,  which  are  highly 
efficient  traps  for  clays  and  suspensions.  Anthropogenic  impact 
using  Pb  isotopes  as  a  tracer  has  so  far  been  documented  for  the 
ancient  harbors  of  Alexandria  (Veron  et  al.,  2006,  2013;  Stanley 
et  al,  2007),  Sidon  (Le  Roux  et  al,  2003),  Marseilles  (Le  Roux 
et  al.,  2005),  and  Rome  (Delile  et  al.,  2014a). 

Applying  similar  methods  to  the  Roman  harbor  of  Ephesus  is 
appealing  because  of  the  status  of  the  Ephesus  city  port  during 
Roman  times  as  an  exceptionally  influential  commercial  and  reli¬ 
gious  center  of  the  ancient  Mediterranean  world.  Ephesus  was  a 
major  town  of  Asia  Minor  and  has  a  long  history  that  began  in  the 
10th  century  BC.  Its  position  at  short  distances  from  both  the  Dar¬ 
danelles  and  the  populated  city  states  of  southern  Greece  gave 
Ephesus  a  strategic  role  in  all  the  wars  affecting  Asia  Minor  and  the 
Aegean  Sea  since  the  Persian  wars  of  the  classical  period.  Its 
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importance  remained  prominent  during  Hellenistic  and  Roman 
times  and  during  the  entire  history  of  the  Byzantine  Empire,  and 
only  declined  as  a  result  of  the  Turkish  conquest.  Because  sedi¬ 
ments  gradually  filled  in  the  inlet  of  the  Cayster  River  (Kii^uk 
Menderes),  the  harbor  of  Ephesus  repeatedly  moved  down  river 
over  the  centuries  (Kraft  et  al,  2000,  2011). 

Here  we  use  samples  from  a  12  m  long  sediment  core  taken  in 
the  Roman  port  of  Ephesus  to  investigate  the  paleo-environmental 
and  hydraulic  evolution  of  the  harbor  using  grain  size  distribution 
analysis,  14C  ages,  major  and  trace  element  geochemistry,  and  Pb 
isotope  compositions.  We  focus  in  particular  on  the  relative  abun¬ 
dances  of  Pb  and  other  chalcophile  elements  in  the  harbor  sedi¬ 
ments  and  discuss  the  respective  status  of  the  anthropogenic  and 
natural  metal  fluxes  and  their  origins  as  deduced  from  the  Pb 
isotope  record. 

2.  Historical  background 

Literature  on  the  history  of  Ephesus  is  abundant  because  of  the 
wealth  of  ruins  left  by  its  different  inhabiting  cultures  and  its  role  in 
the  history  of  this  part  of  the  world  first  as  a  major  religious  center 
dedicated  to  Artemis  and  later  as  one  of  the  leading  churches  of  the 
Mediterranean  world.  For  a  detailed  historical  context  of  the  pre¬ 
sent  work,  the  reader  is  referred  to  the  well-documented  textbook 
by  Foss  (1979)  and  to  Scherrer  (1995).  Here  we  provide  only  a  brief 
overview. 

Different  sites  were  inhabited  in  the  immediate  vicinity  of 
classical  Ephesus  since  the  Neolithic  culture  and  during  the  Bronze 
Age.  The  historical  city  (close  to  the  Artemision)  was  founded  in  the 
10th  century  BC  by  Ionians  and  became  part  of  the  Ionian  League. 
The  classic  site  (at  the  base  of  the  western  side  of  the  Panayirdag) 
was  occupied  around  300  BC  under  Lysimachus,  one  of  Alexander's 
generals,  but  quickly  passed  under  Seleucid  and  then  Ptolemaic 
rules.  After  the  Battle  of  Magnesia  in  190  BC,  Ephesus  came  under 
the  domination  of  Pergamon,  and  finally  became  part  of  the  Roman 
Republic  in  133  BC.  After  the  Mithridatic  wars  (ending  in  63  BC), 
Augustus  made  Ephesus  the  capital  of  Asia  Minor.  At  that  stage,  the 
surface  area  of  the  city,  enclosed  by  the  walls  of  Lysimachus,  is 
thought  to  have  extended  over  more  than  2  km2  and  its  population 
to  have  reached  50,000  inhabitants. 

The  city  and  its  temple  were  destroyed  by  the  Goths  in  262  AD. 
But  Ephesus  was  rebuilt  and  enlarged  by  Constantine  and  soon 
recaptured  most  of  the  importance  it  had  held  since  Hellenistic 
times.  A  burst  of  seismic  activity  between  358  and  365  AD 
repeatedly  destroyed  major  cities  around  the  Aegean  (Guidoboni, 
1994),  including  Ephesus.  In  the  7th  century  AD,  several  addi¬ 
tional  disasters  struck  Ephesus,  notably  the  major  earthquake  of 
614  AD,  as  well  as  the  repeated  sacks  by  Arab,  Frankish,  and  Turkish 
raiders.  Western  Turkey  is  well  known  for  being  subjected  to 
frequent  earthquakes  of  large  magnitude  (e.g.,  Vannucci  et  al., 
2004).  Although  some  dates  are  not  well  established,  particularly 
severe  earthquakes  persistently  ravaged  the  city  in  AD  17,  23,  47, 
178, 194,  262,  275,  337,  358  to  365,  and  614  (Guidoboni,  1994;  Foss, 
1979).  In  AD  1304,  what  was  by  then  left  of  Ephesus  fell  into  the 
hands  of  the  Turks,  and  its  population  was  either  deported  or 
massacred.  These  adverse  troubles  combined  with  the  final  stages 
of  insilting  of  the  harbor  basin,  which  had  incessantly  plagued 
harbor  activity  since  its  early  Hellenistic  days  (Strabo,  XIV.1.24), 
precipitated  the  demise  of  the  harbor  and  the  city  it  served. 

3.  The  study  area 

Ephesus'  harbor  lies  on  the  Aegean  coast  of  Turkey  at  the  western 
extremity  of  the  Kii^uk  Menderes  graben  (KMG)  (Fig.  1).  The  KMG 
corresponds  to  the  catchment  area  of  the  Kii^uk  Menderes  (Cayster) 


river,  which  is  divided  into  five  sub-basins  delimited  by  pre-Miocene 
geology  (Rojay  et  al.,  2005).  The  surrounding  hills  are  composed  of 
crystalline  marble  or  partially  dolomitic  breccias  of  Mesozoic  age 
(Vetters,  1989;  (^akmakoglu,  2007).  The  hills  over  which  Ephesus 
aqueducts  run  also  include  Paleozoic  crystalline  rocks  such  as  gran¬ 
ites,  gneisses,  and  micaschists.  Water  was  brought  to  the  city  by  up  to 
seven  aqueducts  built  between  Archaic  times  and  the  Roman  Empire 
and  repaired  during  different  periods,  notably  after  major  earth¬ 
quakes.  This  point  is  particularly  important  since  all  the  waters  from 
the  aqueducts  terminated  in  the  harbor  where  they  were  susceptible 
to  mixing  with  Cayster  river  water,  marine  water,  and  waste  waters  of 
public  (baths,  fountains)  and  domestic  usage,  as  well  as  with  water 
from  local  workshops  (Ortloff  and  Crouch,  2001 ). 

The  variation  of  relative  sea  level  and  the  westward  migration  of 
the  shoreline  since  Antiquity  have  been  studied  by  Bruckner  (2005) 
and  Pavlopoulos  et  al.  (2012).  Each  century  the  shoreline  pro¬ 
gressed  westward  by,  on  average,  -350  m.  Comparison  of  the 
apparent  sea  level  changes  with  the  values  predicted  by  the 
regional  model  of  Lambeck  and  Purcell  (2005)  indicates  that  sub¬ 
sidence  of  the  coastline  next  to  Ephesus  since  the  classical  period 
was  of  the  order  of  3-7  m. 

Geoarcheological  research  has  been  carried  out  at  the  Ephesus 
site  and  in  the  delta  of  the  Cayster  river  since  the  1990s  (Bruckner, 
1997,  2005;  Kraft  et  al.,  1999,  2000,  2001,  2011;  Stock  et  al.,  2013, 

2014,  unpublished  data).  Besides  reconstruction  of  the  successive 
paleo-environments  and  the  coastline  as  it  has  existed  since 
6000-5000  BC  (Fig.  1),  this  work  also  has  shown  that  delta  pro¬ 
gradation  led  to  multiple  westward  resettlements  of  the  harbor.  The 
ceaseless  fight  against  silting  to  maintain  the  harbor  of  Ephesus  as  a 
functioning  port  during  the  Hellenistic  period  is  first  and  foremost 
reflected  in  the  displacement  of  the  city  to  the  western  side  of  mount 
Pion  (Panayirdag)  by  Lysimachus  in  -290  BC  (Scherrer,  1995). 

4.  Analytical  techniques 

A  sediment  core  about  12  m  long  (EPH  276)  was  drilled  in  the 
hexagonal  Roman  harbor  basin  of  Ephesus  (Fig.  1).  We  sampled  the 
core  at  high  resolution  by  taking  a  total  of  111  samples  (one  sample 
every  10  cm).  The  samples  were  analyzed  for  grain  size  distributions 
(see  Stock  et  al.,  2013,  for  details),  major  and  trace  element  con¬ 
centrations  (Table  SI;  see  Delile  et  al.,  2014b,  for  details),  and  Pb 
isotope  compositions  ( ^able  S2;  see  here  below  and  Delile  et  al., 
2014a,  for  details).  Lead  isotope  compositions  were  obtained  not 
on  the  bulk  sediment,  but,  in  order  to  isolate  potential  anthropo¬ 
genic  components,  on  HBr  leachates.  The  leaching  procedure  con¬ 
sisted  in  first  treating  the  samples  with  chloroform  to  remove  most 
of  the  organic  fraction,  then,  after  rinsing  the  residues  with  clean 
water,  leaching  them  with  dilute  HBr  including  ultrasonicating  and 
heating  steps.  As  shown  for  Portus  (Delile  et  al.,  2014a),  this  tech¬ 
nique  enhances  the  contrast  between  Pb  held  in  surface 
contamination-prone  coatings  and  detrital  silicates.  Carbonates  also 
dissolve  during  the  leaching  process,  but  Pb  contents  of  detrital 
carbonates  are  naturally  low.  As  for  carbonates  precipitated  within 
the  harbor,  they  are  of  biogenic  origin  (cf.  discussion  of  series  D  and  E 
below)  likely  meaning  that  isotope  information  obtained  on 
carbonate-rich  samples  is  consistent  with  that  derived  from  the 
leachates  of  the  rest  of  the  sample  series.  The  HBr  leach  fraction  was 
recovered  for  Pb  separation  by  anion-exchange  chromatography 
using  HBr  as  eluent  of  the  sample  matrix  and  HC1  as  eluent  of  the  Pb. 
Lead  was  also  separated  from  the  residues  of  16  of  the  111  EPH  276 
samples.  The  amounts  of  Pb  extracted  from  all  samples  were  large 
(>1  fig)  and  orders  of  magnitude  above  the  total  procedural  blank  of 
-20  pg.  The  purified  Pb  was  analyzed  for  its  isotopic  composition  by 
multiple-collector  inductively  coupled  plasma  mass  spectrometry 
(Nu  Plasma  500  HR)  at  ENS  Lyon  using  T1  for  instrumental  mass  bias 
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Fig.  1.  Map  of  the  Kii^iik  Menderes  graben  located  on  the  Aegean  coast  of  Turkey  (inset)  with  successive  positions  of  the  shorelines  and  site  of  core  EPH  276  (from  Bruckner,  2005, 
modified  by  Stock  et  al.,  2013  and  this  study). 


correction  and  bracketing  the  samples  with  the  NIST  981  standard 
for  which  the  values  of  Eisele  et  al.  (2003)  were  used. 

Six  samples  were  AMS-14C-dated  ( Table  S3),  complementing  the 
chronostratigraphy  of  neighboring  cores  analyzed  by  Stock  et  al. 
(unpublished  data).  The  Carbon-14  ages  were  obtained  on  frag¬ 
ments  of  wood,  vegetal  matter,  seeds,  and  pollen,  and  are  listed  in 
Table  S3  and  shown  in  Figs.  2  and  3.  Errors  on  raw  radiocarbon  ages 
BP  are  reported  at  the  95  percent  confidence  level  (two  sigma).  The 
measured  14C  (BP)  ages  were  converted  into  BC-AD  dates  relative 
to  the  continental  and  marine  curves  of  Reimer  et  al.  (2009)  using 
the  Clam  software  (Blaauw,  2010). 

Interpretation  of  the  analytical  results  rests  on  different 
methods  of  data  processing.  We  applied  Principal  Component 
Analysis  (PCA)  and  Factor  Analysis  to  major  and  trace  element 
concentration  data,  as  well  as  loss-on-ignition  (L.O.I.)  (Fig.  2).  In  the 
very  large  data  sets  typical  of  those  that  modern  geochemistry  can 
now  produce,  observations  are  often  correlated.  A  common  case  is 
that  of  the  dilution  of  elements  in  sediments  by  detrital  quartz. 
Such  effects  render  the  reading  of  the  underlying  causes  of 
geochemical  variation  and  their  number  difficult.  PCA  consists  in 
rotating  the  data  in  their  multidimensional  space  to  convert  them 
into  uncorrelated  variables  known  as  principal  components.  Un¬ 
correlated  does  not  equate  with  independent,  however,  implying 
that  small  changes  in  rotation  may  affect  all  the  principal  compo¬ 
nents.  PCA  generally  demonstrates  that  the  variability  of  the  ob¬ 
servations  can  be  accounted  for  by  a  very  small  (2-4)  subset  of 
variables  that  carry  the  bulk  of  the  total  variance.  Principal 


components  can  be  calculated  from  the  covariance  matrix  or  from 
the  correlation  matrix.  Factor  Analysis  is  a  related  technique  that 
searches  for  the  minimum  variance  for  an  arbitrary  number  of 
uncorrelated  variables.  It  usually  starts  with  PCA  and  implements 
different  modes  of  rotation  and  weighing. 

In  addition  to  the  PCA  and  Factor  Analysis  we  also  converted  the 
Pb  isotope  compositions  into  their  corresponding  geochemically 
informed  parameters,  which  are  the  model  age  Tmod  and  the 
238u/204pb  (ju)  and  Th/U  (k)  ratios  (Table  S2)  using  the  equations 
given  by  Albarede  et  al.  (2012),  who  also  justified  the  advantages  of 
this  representation  over  those  based  on  raw  Pb  isotope  ratios.  In 
short,  Tmod  is  a  proxy  for  the  tectonic  age  of  crystalline  rocks  and 
their  associated  ore  deposits  (or  depositional  age  for  sediments), 
while  n  is  the  238U/204Pb  and  k  the  Th/U  ratio  of  the  province  in 
which  these  rocks  formed.  Tmod  closely  maps  the  distribution  of  the 
Alpine,  Hercynian,  and  early  Paleozoic  provinces  of  Europe,  while  / x 
delineates  collision  belts,  and  k  is  a  geochemical  parameter  with  a 
remarkable  regional  consistency  related  to  uplift  and  erosion.  Maps 
of  these  parameters  can  be  used  to  divide  Europe  into  coherent 
regions  (Delile  et  al.,  2014a),  which  justify  the  use  of  Tmod,  and  k 
to  determine  provenance  of  archeological  artifacts.  Tmo d,  ft,  and  k  in 
turn  provide  a  rapid  characterization  of  the  geological  environment 
in  which  ores  formed.  A  Matlab  code  is  given  in  Appendix  A  and  an 
Excel  spreadsheet  in  which  to  calculate  these  parameters  will  be 
provided  upon  request.  As  mentioned  above,  rmod  represents  the 
tectonic  age  of  the  geological  province  to  which  a  given  sample 
belongs,  while  jl  is  best  perceived  as  an  indicator  of  whether  this 
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Fig.  2.  Factor  analysis  of  major  and  trace  element  concentrations.  Based  on  Principal  Component  Analysis,  the  number  of  factors  is  limited  to  three.  (A)  Component  loadings.  The 
first  factor  FI  shows  the  trade-off  between  an  (Al,  Ti)-rich  detrital  component  (FI  >  0)  and  a  Ca-rich  carbonate  component  (FI  <  0).  The  second  factor  F2  is  dominated  by  metals  (Ag, 
Pb,  Cu,  Mo)  and  sulfur  and  shows  the  effect  of  anoxia.  The  third  factor  F3  is  dominated  by  the  light  rare-earth  elements  La  and  Ce  and  testifies  to  the  presence  of  heavy  minerals  in 
sand.  (B)  Distribution  of  the  different  factors  with  depth  in  the  column.  The  plots  are  compared  with  the  sedimentary  units  and  the  age-depth  model. 


province  is  a  collision  range  or  a  tectonically  stable  area.  The  vari¬ 
able  k  distinguishes  upper  crust  with  low  k  values  from  middle  and 
lower  crust  with  higher  k  values  (Albarede  et  al.,  2012).  The  pre¬ 
cision  and  accuracy  of  rmod  is  typically  of  a  few  tens  of  Ma,  but, 
occasionally,  the  Tmo(\-jL-K  model  fails  when  the  underlying 
closed-system  assumption  breaks  down  due  to  U  addition  by  recent 
weathering  or  hydrothermal  activity. 

5.  Sedimentary  units  and  the  age-depth  model 

The  core  has  been  divided  into  five  different  units  labeled  A,  B,  C, 
D,  and  E  on  the  basis  of  the  sedimentological  and  geochemical  traits 


described  in  Fig.  3A;  they  span  the  entire  period  of  activity  of  the 
Hellenistic,  Roman,  and  Byzantine  harbor  ( Table  1 ). 

Unit  A  (1200-1080  cm)  exhibits  alternating  brown  and  gray 
varves  composed  of  massive  clayey  silts  with  the  presence  of 
several  beige  to  ocher  fine  layers.  The  C/M  plot  (Fig.  3B)  indicates 
that  the  depositional  processes  are  represented  by  mixed  decan¬ 
tation  and  graded  suspension.  Units  B  and  C  (1080-515  cm)  are 
characterized  by  gray  to  greenish  massive  sandy  silts  with  the 
presence  of  several  beige  to  ocher  fine  layers  enriched  in  sand.  This 
deposit  was  deposited  as  a  graded  suspension  with  embedded  fine 
layers  derived  from  mixed  processes  of  graded  and  uniform  sus¬ 
pension.  Unit  D  (515-290  cm)  consists  of  dark  to  greyish  silts  with 
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Fig.  3.  (A)  Stratigraphy  and  sedimentology  of  core  EPH  276  showing  grain  size  distribution  and  environmental  facies.  (B)  Plot  of  the  grain  size  99  percentile  (D  99)  versus  the 
median  size  (D  50)  for  the  different  samples  analyzed.  The  different  color  groups  correspond  to  different  sedimentation  regimes  (see  legend).  (For  interpretation  of  the  references  to 
color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


variable  clay  (bottom)  and  sand  (top)  enrichments.  From  bottom  to 
top,  deposit  processes  evolve  from  mixed  uniform  and  graded 
suspension  to  a  blend  with  graded  and  uniform  suspensions.  Unit  E 
(<290  cm)  is  composed  of  beige  sandy  silts  with  phragmite  vegetal 
remains.  The  depositional  processes  point  mostly  to  mixed  graded 
and  uniform  suspensions. 

The  age-depth  model  is  based  on  six  14C  ages  (Fig.  4,  Table  S3). 
The  four  oldest  14C  ages  fall  within  a  narrow  time  interval  and  are 


statistically  indistinguishable.  An  approximate  seven  meters  of 
sediment  were  deposited  in  a  few  tens  of  years  during  the  reign  of 
Augustus  or  shortly  after.  Such  an  extraordinarily  fast  sedimentation 
rate  is  consistent  with  a  periodogram  analysis  (e.g.,  Albarede,  1995) 
of  the  magnetic  susceptibility  record.  The  periodogram,  which  is  the 
equivalent  of  a  Fourier  transform  for  unequally  spaced  data,  iden¬ 
tifies  prominent  periodic  fluctuations  in  the  targeted  property,  here 
the  magnetic  susceptibility.  After  removal  of  long-term  variations 
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Table  1 

Depth  range  of  the  sedimentary  units  and  probable  age  assignment. 


Unit 

Depth  (cm) 

Probable  age 

A 

1200-1070 

Roman  Republic 

B 

1070-650 

Early  Roman  Empire 

C 

650-515 

Late  Roman  Empire 

D 

515-300 

Early  Byzantine  (4th— 8th  century) 

E 

300-0 

Late  Byzantine  and  Turkish 

(de-trending)  by  fitting  a  fourth-degree  polynomial,  the  shortest 
values  with  significance  level  P  >  0.95  occur  at  20  and  24  cm  (Fig.  5). 
Longer  wavelength  peaks  probably  reflect  climatic  effects  or  are 
artifacts  of  de-trending.  Assuming  a  seasonal  cause  for  the  observed 
susceptibility  fluctuations  therefore  indicates  a  sedimentation  rate 
of  -20  cm  per  year,  equivalent,  over  the  7  m  of  sediment  with  the 
oldest  14C  ages,  to  -35  years  of  sedimentary  history.  In  contrast,  the 
average  sedimentation  rate  between  the  top  three  14C  samples 
(early  and  late  Byzantine)  is  only  ~0.2  cm  a-1. 

6.  Results  and  discussion 

6.1.  Harbor  hydraulics 

In  order  to  understand  the  hydraulic  dynamics  of  the  harbor,  we 
first  need  to  estimate  its  water  capacity.  The  approximate  di¬ 
mensions  taken  from  aerial  photographs  lead  to  a  volume  of 
(500  x  400  x  5)  m3  =  1.0  x  106  m3  (see  also  Kirbihler,  2013;  Stock 
et  al.,  unpublished  data).  We  assume  that  most  of  the  sedimentary 
layers  were  deposited  during  short  seasonal  flood  events  of  the 
Cayster  River,  whereas  water  running  into  the  harbor  came  from 
different  potential  sources:  seawater,  runoff,  karstic  springs,  and 
aqueducts.  Seawater  and  water  from  springs  and  aqueducts  must 
have  been  largely  clear  of  sedimentary  particles.  Karstic  springs  are 
common  in  the  area  around  the  modern  Lake  Kocagoz  (Somay  et  al., 
2008;  Somay  and  Gemici,  2009),  which  today  exists  on  the  site  of 
the  ancient  harbor  basin.  A  seawater  component  is  present  in  the 
water  from  all  the  lakes,  including  Lake  Kocagoz  (Somay  et  al., 
2008;  Somay  and  Gemici,  2009).  Such  a  component  attests  to  a 
contribution  from  spring  waters  contaminated  by  marine  in¬ 
trusions  into  the  karst.  As  for  other  water  inputs,  Kraft  et  al  (2007) 


Age  BC-AD 

Fig.  4.  Age-depth  model  for  core  EPH  276  deduced  from  the  six  14C  dates  with  ranges 
calculated  using  the  Clam  software  (Blaauw,  2010).  The  size  of  the  data  symbols  re¬ 
flects  the  confidence  level. 


pointed  out  that  all  the  city  sewage  was  diverted  into  the  Great 
Harbor. 

The  up  to  seven  aqueducts  built  during  the  existence  of  the 
Ephesus  port  carried  substantial  volumes  of  water  into  the  harbor. 
Wiplinger  (2013)  quotes  an  estimate  of  0.6  m3  s-1  for  the 
Dergirmendere  aqueduct  alone.  Using  the  model  by  Ortloff  and 
Crouch  (2001),  we  surmise  that  the  total  water  distribution  to  the 
city  from  the  fully-functional  aqueducts  at  the  peak  of  city  prosperity 
may  have  been  over  -2  m3  s-1.  This  number  is  substantial  with 
respect  to  the  mean  discharge  of -11.45  m3  s_1  inferred  for  the  river, 
not  including  flood  events  (Vliegenthart  et  al.,  2007),  indicating  that 
if  other  inputs  such  as  karstic  springs  and  runoff  are  disregarded, 
water  in  the  harbor  was  replaced  by  the  aqueducts  in  merely  six  days. 
This  estimate  is  of  course  an  average  estimate  and  during  seasonal 
droughts  the  ingress  of  seawater  attested  to  by  the  presence  of 
brackish  fauna  also  contributed  to  the  harbor's  overall  water  budget. 

Water  output  is  difficult  to  constrain  independently.  The  harbor 
canal,  whose  construction  may  have  started  as  early  as  during  the 
first  century  AD  when  the  shoreline  swept  past  the  harbor,  was 
narrow  at  the  harbor  entrance  (Kraft  et  al.,  2007).  Assuming  a  cross- 
section  at  the  narrowest  point  of  -50  m2  would  imply  that 
aqueduct-delivered  water  was  leaving  the  canal  at  a  rate  of 
145  m  per  hour,  probably  fast  enough  to  limit  water  ingress  from 
the  sea  under  fair  weather  conditions.  This  velocity  must  have  been 
reduced  by  the  effect  of  draught  and  evaporation  during  the  dry 
(hot)  season,  and  increased  by  local  springs  during  the  wet  (cold) 
season.  The  presence  of  brackish  water  ostracods  and  occasional 
occurrence  of  marine  foraminifera  (Stock  et  al.,  unpublished  data) 
demonstrate  that  the  flow  could  occasionally  be  reversed,  pre¬ 
sumably  as  a  result  of  a  low  water  table  and  reduced  precipitation 
during  the  dry  season.  Over  time,  the  harbor  was  nevertheless 
affected  by  westward  delta  progradation  and  proximity  to  the 
mouth  of  the  Cayster  river  (Fig.  1 ):  by  the  end  of  the  2nd  century  BC, 
the  delta  had  advanced  as  far  as  the  Great  Harbor  (Kraft  et  al.,  2007) 
and  the  canal  had  to  be  extended,  thus  limiting  ingress  of  seawater 
into  the  harbor  basin  even  further. 

To  summarize  harbor  hydrodynamics,  the  ‘normal’  situation  is 
that  of  a  basin  steadily  filled  by  polluted  urban  water  initially  brought 
to  the  harbor  by  aqueducts  and  local  springs  and  quickly  evacuated 
through  a  canal  with  little  ingress  from  the  sea.  As  long  as  the 
coastline  is  not  too  distant,  some  seawater  may  be  admitted  during 
the  dry  season,  while  floods  of,  in  the  present  case,  the  Cayster  River 
dominated  the  water  balance  during  spells  of  heavy  rain.  Silting  in  of 
the  harbor  would  have  been  caused  only  by  floods,  which  today  are 
known  to  carry  up  to  100-150  m3  s_1  of  water  (Vliegenthart  et  al., 
2007).  The  Romans  went  to  great  length  to  protect  the  harbor  from 
river  floods.  Kraft  et  al.  (2000,  2011)  mention  that,  in  the  early  2nd 
century  AD,  Hadrian  sought  to  divert  the  Cayster  River  with  an  18  m 
high  dam  and  also  made  multiple  attempts  to  dredge  the  harbor. 

6.2.  Environmental  conditions  in  the  harbor  basin 

In  order  to  assess  the  environmental  conditions  that  prevailed 
during  sedimentation,  we  plotted  the  concentrations  of  first-row 
transition  elements  (Ti-Zn)  and  other  metallic  elements  (Ga,  Pb, 
Mo,  Bi,  Cd,  Ag,  As,  and  Sb)  normalized  to  the  upper-crust  concen¬ 
trations  of  Rudnick  and  Gao  (2003)  (Fig.  SI).  Factor  analysis  of 
major  and  trace  element  abundances  leads  to  the  identification  of 
three  major  components. 

1.  The  first  factor  opposes  elements  indicative  of  the  detrital  load  of 

the  river  (Al,  Ti,  Mg,  etc.)  to  those  distinctive  of  carbonate 

minerals  (Ca,  Sr)  and  L.O.I. 

2.  The  second  factor  singles  out  chalcophile  elements  that,  as 

attested  to  by  the  presence  of  sulfur  in  this  group,  precipitate  as 
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Fig.  5.  Periodogram  of  magnetic  susceptibility  in  core  EPH  276  between  440  and  1100  cm.  Mean  sampling  interval  is  generally  ~1  cm.  Four  extreme  values  corresponding  to 
discontinuities  were  removed  from  the  data  set.  A  fourth-order  polynomial  was  then  fitted  to  the  data  to  remove  the  long-term  trend.  Peaks  correspond  to  dominant  periods,  with 
confidence  levels  P  in  percent  in  parentheses.  We  consider  that  the  shortest  periods  with  P  >  95%  correspond  to  dominant  annual  varves,  while  longer  periods  correspond  to 
climatic  effects  or  to  de-trending  artifacts.  The  periodogram  is  interpreted  as  indicating  an  average  sedimentation  rate  of  20  cm  per  year. 


sulfide  under  anoxic  conditions  (Pb,  Ag,  Cu,  Ni,  Mo),  or  are 
particularly  sensitive  to  redox  conditions  (U,  Cr).  When  the  el¬ 
ements  embedded  in  this  factor  are  normalized  to  Al  (Fig.  6),  as  a 
means  of  accounting  for  the  variable  abundance  of  the  detrital 
component,  and  plotted  against  depth  in  the  core,  a  sharp  in¬ 
crease  is  observed  at  520  cm  depth.  The  significance  of  this 
factor  deserves  some  discussion  because  Pb,  Ag,  and  Cu  may  also 
be  seen  as  representing  an  anthropogenic  component.  Fig.  6  and 
SI  further  show  the  striking  consistency  of  these  metals  among 
themselves  and  with  respect  to  sulfur.  Such  regular  behavior  is 
not  supportive  of  random  contamination  by  a  particular  metal, 
such  as  Pb.  The  Mo-Pb  correlation  is  very  strong  (r  =  0.90 
excluding  the  top  five  samples  likely  contaminated  by  gasoline 
Pb)  as  is  the  Ag-Pb  correlation  (r  =  0.95).  This  factor  therefore 
reflects  more  on  changing  redox  conditions  in  the  harbor  than 
on  anthropogenic  pollution. 


3.  The  third  factor  is  dominated  by  La,  Ce  and,  to  a  lesser  extent,  Se. 
Most  other  loadings  are  very  small,  except  possibly  P.  The  weak 
negative  correlation  between  excess  La  and  Ce  on  the  one  hand 
and  P  deficit  on  the  other  hand  suggests  the  presence  of  non- 
phosphatic  rare-earth  minerals,  such  as  allanite,  notably  in  the 
coarse  silts  between  515  and  650  cm. 

The  accumulation  of  so  much  sediment  in  a  matter  of  decades 
requires  an  explanation,  especially  since  the  thickness  of  the  newly 
deposited  layers  exceeds  the  water  depth  usually  assumed  for  the 
harbor  (4-6  m),  even  next  to  the  mole  (Stock  et  al.,  unpublished 
data).  One  factor  clearly  is  the  westward  progradation  of  the 
shoreline  past  the  harbor  at  about  the  time  of  fast  sedimentation.  F. 
Stock  (personal  communication)  obtained  a  14C  age  of  44  BC-AD  52 
for  the  silting  in  of  the  harbor  canal  consistent  with  the  present 
finding.  The  sharp  geochemical  discontinuities  at  650  and  550  m 
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Fig.  6.  Downcore  variations  of  calcium  and  metal  concentrations  normalized  to  aluminum.  White-  and  gray-shaded  bands  delineate  the  stratigraphic  units  A— E,  with  unit  A  being 
deepest  and  unit  E  shallowest.  Note  the  discontinuities  at  -1080  and  -550  cm  depth,  notably  the  increase  in  S  and  chalcophile  elements  at  the  latter. 
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are  flagged  by  strong  peaks  of  magnetic  susceptibility  (Fig.  3A).  The 
efficiency  of  sediment  confinement  by  the  harbor  prior  to  the  three 
lowermost  units  is  staggering,  while  the  sudden  drop  in  sedimen¬ 
tation  rate  and  the  short  but  intense  episodes  of  high  magnetic 
susceptibility  require  the  intervention  of  a  brutal  event.  A  probable 
cause  for  this  discontinuity  is  an  abrupt  jump  of  the  Cayster  River 
channel  triggered  by  the  abandonment  of  a  meander  (avulsions)  or 
by  exceptional  floods  (Brown,  1997).  Co-seismic  vertical  move¬ 
ments  (Pavlopoulos  et  al.,  2012)  associated  with  the  major  AD  17, 
23,  and  47  earthquakes  may  also  have  played  a  role. 

The  lowermost  unit  A  (1200-1080  cm)  was  deposited  during  the 
Roman  Republic.  It  is  consistently  dominated  by  silt  (FI  >  0)  with 
anoxic  influence  (F2  <  0).  The  anoxic  conditions  of  the  basin  bottom 
as  attested  to  by  abundant  S,  Mo,  and  U  (Fig.  2),  small  excesses  of 
Mo  and  Ag  (Fig.  SI),  and  persistence  of  seasonal  varves,  indicate 
that  the  terrigenous  flux  into  the  early  Roman  harbor  of  Ephesus 
during  the  1st  century  BC  was  not  noticeably  oxidized  whether 
water  input  was  freshwater  or  seawater.  Grain  size  distribution 
(Fig.  3)  reflects  an  environment  where  decantation  is  important 
(Bravard  and  Peiry,  1999;  Bravard  et  al.,  2014).  Input  of  oxygenated 
freshwater  into  the  harbor,  regardless  of  its  source,  therefore  was 
limited  and  whatever  water  was  added  by  the  aqueduct  system 
must  have  been  dominated  by  sewage. 

Unit  B  (1080-650  cm)  continues  to  show  the  prevalence  of  the 
detrital  flux  (FI  >  0),  but  now  with  evidence  of  oxygenation 
(F2  >  0).  The  transition-element  pattern  typically  is  crustal  in  origin 
and  no  visible  anomaly  of  Mo  and  Ag  is  observed  (Fig.  SI).  Grain- 
size  analysis  indicates  graded  grain  size  distributions  by  turbulent 
waters,  reflecting  that,  even  at  times  of  flood,  water  was  being 
constantly  evacuated  from  the  harbor.  Ephesus  counted  up  to  seven 
aqueducts  implying  that  the  early  Roman  harbor  was  saved  from 
silting  in  as  much  by  water  from  its  many  aqueducts  continuously 
flushing  the  basin  as  by  the  Roman  engineers.  As  shown  by  the 
return  of  some  decantation  events  (Fig.  3),  the  aqueducts  made 
silting  depend  on  a  fully  functioning  water  distribution  system.  In 
this  respect,  the  Menderes  area  is  seismically  active  ( Vannucci  et  al., 
2004)  and  major  earthquakes  were  particularly  disruptive  to  the 
long  and  complex  Ephesus  aqueduct  network  (Passchier  et  al., 


2013).  Reduction  of  the  water  input  by  the  seismic  destruction  of 
aqueducts  translates  into  reduced  water  egress  from  the  harbor 
basin  and  hence  enhanced  efficiency  of  its  role  as  a  sediment  trap. 
Silting  in  of  the  harbor  in  the  aftermath  of  major  earthquakes 
therefore  became  collateral  damage  to  the  rest  of  the  disasters 
caused  by  the  seismic  activity. 

Transition  to  unit  C  (Fig.  2)  (650-550  cm)  is  heralded  by  a 
peak  of  magnetic  susceptibility  (Fig.  3A).  Highly  negative  values 
of  factor  3,  i.e.,  higher  Se,  La,  and  Ce  contents,  reflect  lesser 
dilution  of  minor  elements  by  quartz  and  carbonate.  The  varia¬ 
tion  patterns  of  transition  elements  and  other  metals  are  very 
similar  to  those  of  unit  B.  As  already  observed  for  rivers  (Yang 
et  al.,  2002),  a  strong  correlation  exists  between  grain  size  and 
lanthanide  concentrations  (Zhang  et  al.,  1998;  Yang  et  al.,  2002). 
This  geochemical  change  is  consistent  with  a  sand  fraction  in  unit 
C  smaller  than  that  in  unit  B.  Unit  C  shows  some  transient 
geochemical  features  (Fig.  7),  true  harbingers  of  the  major 
changes  that  would  profoundly  affect  unit  D,  notably  an  increase 
in  sulfur  and  heavy  metal  contents. 

The  transition  (550-515  cm)  between  unit  C  and  unit  D 
(515-300  cm)  also  is  announced  by  a  strong  peak  of  magnetic 
susceptibility  (Fig.  3A),  corresponding  to  a  strong  compositional 
shift  with  a  surge  of  the  biogenic  component  (FI  <  0)  due  to 
degraded  ventilation  of  bottom  waters  by  eutrophication  (F2  <  0). 
Sulfide  reduction  and  precipitation  is  attested  to  by  a  sudden  two- 
order-of-magnitude  increase  in  the  S/Al  ratio  (Fig.  6).  The  surge  in 
sulfur,  Zn,  Ni,  and  Co  conspicuously  follows  the  surge  of  Ca,  Pb,  Ag, 
Cu,  Cd,  Mo,  and  Cr  by  some  30  cm  in  the  core.  This  delay,  which  may 
have  been  as  short  as  a  few  years  and  possibly  was  only  one  or  two 
years,  is  visible  in  the  plot  of  Fig.  7  as  a  pronounced  negative 
excursion  of  ratios  such  as  Ca/S,  Mo/S,  and  Pb/S  between  520  and 
550  cm  depth.  These  characteristics  together  with  high  Sr  abun¬ 
dance  and  the  presence  of  fine  calcareous  layers  (Kylander  et  al., 
2011;  Martin-Puertas  et  al.,  2011)  show  that  sulfide  precipitation 
predated  the  development  of  eutrophic  conditions  manifested  by 
the  rise  in  Ca  and  was  due  to  the  sudden  isolation  of  the  harbor 
from  ventilated  waters.  The  trend  of  decreasing  ratios  of  chalco- 
phile  elements  relative  to  sulfur,  which  was  perceptible  already  in 


100 


i/> 

"?u 

+■* 

0) 


0.01 


0.001 


0.0001 


0.00001 


0.000001 


—  Zn/S 

—  Ni/S 

—  Co/S 

—  Ca/S 

—  Pb/S 

—  Ag/S 

—  Cu/S 

—  Cd/S 


100  200  300  400  500  600  700  800  900  1000  1100  1200 


Depth  below  surface  (cm) 


Fig.  7.  Downcore  variations  of  calcium  and  metal  concentrations  normalized  to  sulfur  (see  caption  of  Fig.  4  for  details).  Two  major  discontinuities  are  again  observed  at  -1080  and 
-550  cm  depth,  which  attest  to  sudden  changes  in  sediment  oxygenation.  The  1080  cm  discontinuity  reflects  the  improved  oxygenation  of  the  harbor  in  the  early  1st  century  AD. 
Most  elements  show  a  negative  excursion  between  550  and  520  cm  indicating  that  the  rise  in  sulfur  precedes  the  rise  in  Ca  and  most  other  metals;  exceptions  are  Zn  and  Co,  which 
are  in  phase  with  S.  Although  a  slow  trend  towards  anoxia  can  be  seen  in  sediments  below  the  550  cm  level,  the  rate  of  ventilation  is  greatly  reduced  after  this  episode. 
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Fig.  8.  Downcore  variations  of  Pb/Alnorm,  Tmod,  g  (238U/204Pb),  and  k  (Th/U)  compared  with  variations  in  Ca/Al  and  Mo/Al.  Pb/AI  is  normalized  to  the  upper  continental  crust  average 
of  McLennan  (2001).  Major  rmod  discontinuities  are  observed  between  units  C  and  D,  and  at  the  base  of  unit  E.  A  major  discontinuity  in  both  g  and  k  is  observed  between  units  A  and 
B.  The  variability  within  unit  E  reflects  an  anthropogenic  component  of  uncertain  origin. 


the  early  harbor,  markedly  changes  at  the  C-D  boundary,  and  the 
rate  at  which  these  ratios  change  significantly  increases  as  well. 
Again,  the  correlation  between  Cr,  Cu,  Pb,  Mo,  and  Ag  (Figs.  6  and  7, 
SI)  excesses  are  not  in  favor  of  selective  pollution  by  metallurgical 
or  any  other  industrial  activities.  The  destruction  of  the  aqueducts 
by  the  major  earthquakes  ravaging  the  city,  such  as  the  AD  17,  23, 
and  47  events  (see  a  discussion  in  Guidoboni,  1994),  and  the  AD  64 
AD  cleaning  credited  to  Barea  Soranus  by  Tacitus  (XVI, 23)  may  have 
combined  with  the  increasing  silting  of  the  harbor  entrance  upon 
westward  progradation  of  the  delta  (Fig.  1)  (Bruckner,  2005)  to 
modify  the  hydraulic  regime  of  the  harbor.  The  top  of  unit  D  records 
a  short-lived  return  of  better  oxygenated  conditions  which,  with 
the  caution  due  the  age-depth  model,  may  correspond  to  the 
revival  of  the  harbor  by  Justinian  (early  and  mid-  6th  century;  Foss, 
1979;  Scherrer,  1995). 

The  age-depth  model  (Fig.  4)  places  the  transition  between 
units  D  and  E  (-300  cm)  in  the  9th  century.  Carbonate  precipita¬ 
tion  dominates  unit  E  as  it  did  the  lower  part  of  unit  D  (FI  <  0) 
indicating  a  negative  water  balance  (Fig.  2B;  Martin-Puertas  et  al., 
2011;  Delile  et  al.,  2014b).  Excesses  of  Cr,  Cu,  Mo,  and  Ag  are  still 
well  correlated  (Fig.  2A)  and  confirm  the  persistence  of  a  sulfur- 
rich,  oxygen-deficient  eutrophic  regime,  but,  as  shown  by  the 
positive  F2  values,  with  oxygen  deficiency  being  less  pronounced 
than  in  the  underlying  unit  D.  The  water  deficit  caused  conditions 
to  evolve  towards  a  peatland  environment  consistent  with  the 
considerable  extension  of  the  Cayster  delta  at  this  time  (Fig.  1). 
The  modern  estuary  of  the  Kii^tik  Menderes  is  wetland  dotted 
with  alkaline  lakes  recharged  from  precipitation  and  local  karstic 
springs.  In  late  Byzantine  times,  the  harbor  may  have  been  func¬ 
tional,  but  appears  to  have  been  increasingly  cut  off  from  the  sea 
and  the  river  (Kraft  et  al.,  2011)  (Fig.  1).  Some  of  the  shallow  core 


samples  show  excess  Pb  of  probable  but  uncertain  anthropic 
origin  (Fig.  8). 

To  sum  up  on  environmental  conditions,  core  EPH-276  holds  the 
record  of  anoxic  conditions  prevailing  at  times  in  the  harbor,  likely 
compromising  the  control  of  harbor  hydrodynamics  by  human 
activities.  Lead  isotopes  are  expected  to  shed  light  on  the  magni¬ 
tude  of  anthropogenic  contamination  at  the  time  of  sediment 
deposition  and  this  is  what  the  next  section  will  be  addressing. 

6.3.  Interpretation  ofPb  abundance  and  isotopic  signals 

Here  we  focus  on  the  Pb  isotope  compositions  of  the  leached 
fractions  only  because  this  is  where  chances  of  observing  anthro¬ 
pogenic  input  are  maximum.  Fig.  8  shows  the  Pb/AI  ratio  together 
with  the  Pb  isotope  data  in  the  form  of  three  geochemical  parame¬ 
ters,  Tmod  (Ma  =  million  years),  p,  and  k.  Based  on  these  parameters 
and  the  enrichment  factor  of  Pb  as  represented  by  Pb/AI  and  closely 
tracking  F2,  the  chronostratigraphic  evolution  of  these  four  curves 
shows  remarkable  discontinuities  (Fig.  2).  In  agreement  with  what 
was  discussed  above  for  other  metals,  the  transition  between  units  A 
and  B  stands  out  clearly  by  a  marked  drop  in  Pb/AI.  A  subtle  increase 
in  g  and  k ,  while  Tmod  remains  young  (~80  Ma  or  Upper  Cretaceous), 
is  evidence  of  change  in  the  sources  of  Pb.  The  next  discontinuity 
takes  place  between  units  C  and  D.  The  increase  in  Pb/AI  and  Ca/Al  is 
associated  with  older  Tmod  (-120  Ma  or  Early  Cretaceous)  and  lower 
values  of  g  and  k.  The  Pb/AI  ratio  decreases  steadily  throughout  unit 
D  regardless  of  the  changes  in  the  major  Ca/Al  dip  at  3.8  m  depth  that 
we  assigned  above  to  the  6th  century.  From  unit  D  to  E,  most  Tmod 
values  exceed  240  Ma  and  the  k  values  decrease  below  the  level  of 
previous  values.  The  samples  at  the  top  of  the  core  seem  to  be  largely 
influenced  by  a  modern  anthropogenic  component. 
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The  hydraulics  of  the  harbor,  notably  its  volume  and  output, 
may  affect  harbor  oxygenation  and  thereby  the  metal  contents  of 
sediments.  In  contrast,  changing  Pb  isotope  compositions  require 
changes  in  the  relative  contribution  of  all  the  sources  of  this  metal. 
Lead  isotopes  reveal  the  nature  and  relative  strengths  of  the 
following  potential  sources: 

1.  Local  natural  sources ,  which  are  multiple.  Lead  from  the  Cayster 
River  comes  during  short-lived  seasonal  flood  events.  Seawater 
should  be  extremely  poor  in  Pb,  while  brackish  water  from  the 
estuary  should  be  Pb-depleted  by  iron  flocculation  in  the  mixing 
zone.  A  contribution  from  the  runoff  and  from  karstic  springs 
that  discharge  from  marble-schists  and  marble-alluvium  con¬ 
tacts  (Somay  et  al.,  2008;  Somay  and  Gemici,  2009)  should  also 
be  considered. 

2.  Lead  from  the  main  water  distribution  system.  Such  a  component 
may  come  from  the  underground  of  the  springs.  It  can  also  be  ac¬ 
quired  during  transit  from  the  aqueduct  masonry,  which  includes 
mortar  produced  from  local  limestones.  The  laminated  deposits 
observed  in  some  aqueducts  (Passchier  et  al.,  2011,  2013)  indicate 
hard  water  with  excess  alkalinity,  which  does  not  favor  the  idea 
that  Pb  was  leached  out  of  these  conduits.  The  seemingly  high 
concentrations  of  chalcophile  elements  in  the  sediments,  first  and 
foremost  Pb  but  also  Co  and  Mo,  do  not  entail  a  pollution  signal. 

3.  The  network  of  aqueducts ,  which  is  unlikely  to  be  a  major  Pb 
contributor.  The  secondary  water  distribution  system  at  Ephe¬ 
sus  is  dominated  by  terra  cotta  pipes  and  only  rarely  involves 
small-diameter  lead  pipes  or  fistulae  (Ortloff  and  Crouch,  2001 ). 

4.  Anthropogenic  Pb  from  local  workshops  or  ballast  dropped  by 
merchant  ships.  Zabehlicky  (1995)  writes  that  a  lead  anchor 
14.2  cm  long  was  found  during  excavations,  as  well  as  lead 
interpreted  as  ballast,  which  as  much  as  hinting  at  a  potential 
source  of  pollution,  signals  that  the  dissolution  of  Pb  artifacts  is  an 
exceedingly  slow  process.  The  presence  of  an  arsenal  on  the 
harbor  site  at  the  time  of  Augustus  was  noted  by  Strabo  (XIV.1.24). 
There  is  no  doubt  that  a  city  with  the  population  of  Ephesus  at  its 
best  periods  would  have  to  rely  on  local  metallurgy.  In  sediments, 
however,  a  geochemical  signal  of  pollution  is  difficult  to  detect, 
and  to  which  extent  Pb  artifacts  attest  to  wholesale  contamina¬ 
tion  of  harbor  sediments  is  not  clear.  The  stability  of  Pb  isotope 
compositions  over  long  periods  of  time  (decades,  even  centuries) 
does  not  bring  to  mind  metal  supply  in  troubled  times.  Even  more 
conclusive  is  the  observation  that,  except  for  unit  E,  Cr,  Cu,  Mo,  Ag, 
and  Cu  coherently  track  Pb  and  Ca:  although  these  elements  are 
sensitive  to  redox  conditions,  they  were  not  involved  in  Pb  met¬ 
allurgy.  The  Pb/Ag  ratio  remains  remarkably  stable,  while  the 
record  of  Ni  and  Co,  for  which  extractive  metallurgy  was  un¬ 
known  at  the  time  and  which  additionally  are  not  found  in  the 
same  ores  as  Pb,  also  follow  the  Pb  record  with  depth.  The  lesser 
impact  of  anthropogenic  Pb  pollution  at  Ephesus  relative  to  Por- 
tus  (Delile  et  al.,  2014a)  directly  reflects  that  the  urban  water 
distribution  systems  used  different  materials,  terra  cotta  for  the 
former  (Ortloff  and  Crouch,  2001)  and  lead  fistulae  for  the  latter. 

The  Cretaceous  model  ages  of  80  Ma  and  120  Ma  observed  in  the 
lower  part  of  the  core  are  consistent  with  those  of  the  carbonate 
hills  surrounding  Ephesus  and  may  simply  register  Pb  from  the 
water  distribution  system  (natural  springs  or  conduits).  This 
interpretation  is  supported  by  the  similar  Pb  isotope  compositions 
of  leachate-residue  pairs  (  fable  S2).  The  residues  of  the  samples 
with  Tmo d  values  >200  Ma  have  not  been  analyzed,  but  these  ages 
are  consistent  with  the  age  of  the  Paleozoic  and  Triassic  crystalline 
basement  of  the  Menderes  Massif  (Vetters,  1989;  Bozkurt,  2007; 
(Jakmakoglu,  2007;  Giirer  et  al.,  2009)  (see  maps  in  the  supple¬ 
mentary  material  of  Delile  et  al.,  2014a). 


The  appearance  of  old  Tmod  values  shortly  after  the  beginning  of 
unit  E,  which  have  not  been  observed  at  earlier  times,  coincides 
with  the  onset  of  carbonate  sedimentation  (up  to  45  wt.%  CaO). 
Runoff  and  karstic  springs  more  or  less  contaminated  by  seawater 
(Somay  et  al.,  2008;  Somay  and  Gemici,  2009)  are  therefore  left  as 
the  main  steady  sources  of  water  in  the  ancient  port,  which  since 
has  become  the  modern  Lake  Kocagoz.  Two  competing  in¬ 
terpretations  are  left:  (1 )  while  ancient  aqueducts  were  bringing  in 
Cretaceous  Pb,  runoff  and  karstic  springs  now  bring  in  Paleozoic  Pb; 
or  (2)  old  Tmod  values  reflect  some  anthropic  influence  of  poorly 
constrained  origin.  The  younger  samples  with  Tmod  values  >200  Ma 
may  represent  sources  in  the  Menderes  region,  but  Pb  from  Thrace, 
or  even  from  western  Europe,  cannot  be  excluded. 

Historical  evidence  in  favor  of  Ephesus  hosting  significant  in¬ 
dustrial  activity  in  the  11th  and  12th  centuries  (Foss,  1979,  pp. 
120-123),  however,  is  faint.  Foss  (1979,  p.  113)  further  argues  that 
from  the  8th  century  onwards  the  harbor  district  was  literally 
abandoned.  The  center  of  town  moved  to  the  hill  of  Ayasuluk,  while 
new  ports,  such  as  Scalanova  (ancient  Phygela)  built  on  the  site  of 
modern  Ku^adasi,  gradually  took  over  the  silted  in  harbor  of  Ephesus. 

Nevertheless,  from  the  middle  of  the  9th  to  the  10th  century,  the 
victories  of  the  Byzantines  against  the  Arabs  in  Asia  Minor  enabled 
Ephesus  to  regain  a  preeminent  position  in  the  Empire  (Foss,  1979). 
A  phase  of  economic  development  would  be  consistent  with  the 
influx  of  Hercynian  Pb  into  the  most  recent  harbor  deposits.  The 
10th  century  medieval  economic  revolution  in  Europe  favored  the 
development  of  trade  between  Europe  and  the  Orient.  West  or 
north  European  sources  of  Pb  cannot  be  excluded  for  that  period. 

7.  Concluding  remarks 

Major  changes  in  the  lithology,  grain-size  distribution,  major 
and  trace  element  chemistry,  and  Pb  isotope  compositions  of  the 
harbor  sediments  at  Ephesus  reflect  the  history  of  the  water  dis¬ 
tribution  system  of  this  port,  notably  in  response  to  the  increasing 
and  declining  needs  of  a  population  inhabiting  a  city  that  at  several 
points  in  history  was  one  of  the  largest  of  the  Eastern  Mediterra¬ 
nean.  Throughout  its  history,  the  Ephesus  port  was  affected  by 
major  disruptive  events  in  the  form  of  earthquakes  and  invasions, 
both  of  which  were  particularly  effective  at  destroying  aqueducts. 

Progressive  silting  in  of  the  harbor  responded  to  the  westward 
migration  of  the  coastline  and  to  human  maintenance  aimed  at 
keeping  the  harbor  functional.  A  single  major  disruptive  event  located 
at  550  cm  core  depth  and  heralding  a  two-order-of-magnitude  drop 
in  sedimentation  rate  and  the  development  of  anoxia  in  the  harbor  is 
clearly  visible  in  the  major  and  trace  element  record.  Although  this 
event  may  have  unidentified  military  or  seismic  causes,  we  favor  a 
durable  displacement  of  the  river  course,  which  starved  the  harbor 
from  further  silt  input.  Overall,  despite  the  presence  of  metallic  arti¬ 
facts  in  the  harbor,  the  record  of  metal  concentrations,  in  particular 
the  Pb  isotope  record,  suggests  that  pollution  of  the  harbor  was 
subdued  relative  to  other  inputs,  notably  those  of  aqueducts,  except 
near  the  time  of  harbor  abandonment  (unit  E). 

Dating  and  identifying  the  seemingly  key  event  located  in  the 
present  sediment  core  at  550  cm  depth,  as  well  as  in  other  cores 
from  the  same  basin,  is  a  new  and  major  challenge.  This  event 
conspicuously  marks  the  end  of  the  dynamic  regime  controlling  the 
harbor  water  during  all  of  the  Roman  Empire  and  clearly  represents 
a  major  disturbance  in  the  history  of  life  in  Ephesus. 
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Appendix  A.  Matlab  code  for  Tmo d»  and  k  calculations 


The  input  data  should  be  provided  as  an  Excel  file  input.xlsx  made  of  three  columns, holding  1()6Pb/204Pb, 

207Pb/204Pb,  and  208Pb/204Pb. 


%  Matlab  code  for  model  age  calculations  (Francis  Albarede) 

% 

A=xlsread('input.xlsx',  1 ); 

m=size(A,l);  %  find  m  the  number  of  samples 


%  decay  constants 

18=0. 1 55 125;15=0.98485;12=0. 049475; 

%  common  Pb  parameters  from  Albarede  and  Juteau  (1984) 
xstar=l 8.7500;ystar=l 5.63;zstar=38.83;  kappastar=3.90;mustar=9.66; 

%  Initialize 

FF=zeros(m,  1  );mu=zeros(m,  1  );dmu=zeros(m,  1  );kappa=zeros(m,  1  );dkappa=zeros(m,  1  );T  mod=zeros(m,  1  );Tinit= 

i; 

T0=3.8;  %  Beginning  stage  2  (Ga),  see  Stacey  and  Kramers  (1975) 
options=optimset(’DisplayVoff); 
for  i=l  :m  %  loop  through  all  samples 
x=A(i,l);  y=A(i,2);  z=A(i,3); 
myf=@(T)TmPb(T,TO,x,y); 

[Tmod(i),FF(i),exitflag(i,l)]=fsolve(myf,Tinit, options);  %  use  Matlab  fsolve  with  function  TmPb 
Tl=Tmod(i); 

dmu(i)=(x-xstar+mustar*(exp(18*Tl)-l))/(exp(18*T0)-exp(18*Tl)); 
mu(  i  )=mustar+dmu(  i ) ; 

dkappa(i)=(z-zstar+mustar*kappastar*(exp(12*T  1 )- 1  )-kappastar*dmu(i)*(exp(12*T0)- 
exp(12*T  1  )))/(exp(12*T0)-exp(12*T  1  ))/mu(i); 

kappa(i)=kappastar+dkappa(i); 

end 

B=[A,  1 000*Tmod, mu, kappa, FF,exitflag] 

%  B  holds  columnwise  the  original  ratios  (A),  the  model  age  (in  Ma),  mu,  kappa,  the 
%  exit  value  of  the  function  to  solve  (should  be  less  than  le-9)  and  an 
%  exit  parameter  ( 1  expected) 

function  F  =  TmPb(T,TO,x,y) 
mustar=9.66; 

xstar=l  8.7500;ystar=15.63;  %AJ 
18=0.1 55 125;15=0. 98485; 

p0=(exp(15*T0)-exp(15*T))/(exp(18*T0)-exp(18*T))/137.8; 
p  l=(exp(15*T)- 1  )/(exp(18*T)- 1  )/l  37.8; 

% 

if  x~=xstar  %  Eq.  12  in  Albarede  et  al.  (2012) 

F=(y-ystar)/(x-xstar)-p0-mustar*(exp(18  *T)- 1  )*(p0-p  1  )/(x-xstar); 
else 

F=(y-ystar)-mustar*(exp(18*T)-l  )*(p0-pl ); 
end 
end 
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Appendix.  BSupplementary  data 

Supplementary  data  related  to  this  article  can  be  found  at  http:// 
dx.doi.org/10.1016/j.jas.2014.10.002. 
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